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Wc study the processes e+e^ 3(7r+7r~)7, 2(7r+7r^ 77*^)7 and K'^ K~ 2(tx^ -k')") , with the photon 
radiated from the initial state. About 20,000, 33,000 and 4,000 fully reconstructed events, respec- 
tively, have been selected from 232 fb^^ of BABAR data. The invariant mass of the liadronic final 
state defines the effective e'^e" center-of-mass energy, so that these data can be compared with 
the corresponding direct e+e" measurements. Prom the 3(7r"'"7r~), 2(7r"'"7r~7r°) and K~2{-k^-k~) 
mass spectra, the cross sections for the processes e"'"e~ — > 3(7r"'"7r~), e"'"e~ — » 2(7r"'"7r~7r°) and 
e'^e~ — > K'^ K~ 2{tt'^tx~) are measured for center-of-mass energies from production threshold to 
4.5 GeV. The uncertainty in the cross section measurement is typically 6-15%. We observe the J/ij) 
in all these final states and measure the corresponding branching fractions. 
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PACS numbers: 13.66.Bc, 14.40.Cs, 13.25.Gv, 13.25.Jx, 13.20.Jf 



I. INTRODUCTION 

The idea of utilizing initial-state radiation (ISR) from a 
high-mass state to explore electron-positron processes at 
all energies below that state was outlined in Ref. [1] . The 
possibility of exploiting such processes in high luminosity 
(j)- and S-factories was discussed in Refs. [2-4] and moti- 
vates the hadronic cross section measurement described 
in this paper. This is of particular interest because of the 
small deviation of the measured muon g — 2 value from 
that predicted by the Standard Model [5], where hadronic 
loop contributions are obtained from e+e" experiments 
at low center-of-mass (cm.) energies. The study of ISR 
events at i3-factories provides independent and contigu- 
ous measurements of hadronic cross sections in this en- 
ergy region and also contributes to the investigation of 
low-mass resonance spectroscopy. 

The ISR cross section for a particular hadronic final 
state / is related to the corresponding e+e~ cross section 
CT/(s) by: 



daf{s, x) 
dx 



= W{s,x)-af{s{l-x)) , 



(1) 



where x = 2Ey/^; is the energy of the ISR photon in 
the nominal e'^e~ cm. frame; -^i = Ec.m. is the nominal 
e+e~ cm. energy; and i/s(l — x) is the effective cm. en- 
ergy at which the final state / is produced. The invariant 
mass of the hadronic final state is used to measure the 
effective e+e^ cm. energy. The function W{s,x) is cal- 
culated with better than 1% accuracy (see for example 
Ref. [4]) and describes the probability density function 
for ISR photon emission. ISR photons are produced at 
all angles, with a distribution peaking at small angles 
with respect to the axis of the beams, and are required 
to be detected in the electromagnetic calorimeter (EMC) 
of the BABAR detector. The acceptance for such photons 
is 10-15% [4] depending on applied selections. 

An important advantage of ISR data is that the entire 
range of effective cm. energies is scanned in one experi- 
ment. This avoids the relative normalization uncertain- 
tics that inevitably arise when data from different exper- 
iments, or from different machine settings, are combined. 

A disadvantage of the ISR measurement is that the 
mass resolution is much poorer than can be obtained in 
direct annihilation. The resolution and absolute energy 
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scale can be monitored directly using the measured width 
and mass of the J/ip resonance produced in the reaction 
e'^e~ J/ip"f- Using a kinematic fit to this reaction, we 
find the resolution to be about 9 MeV/c^ for decays of 
J/t/j in the 3(7r+7r^) mode and about 15 MeV/c^ in the 
2(7r+7r^7r") mode as will be shown later. 

Studies of e+e^ ^ /i+/i^7 and several multi-hadron 
ISR processes using BABAR data have been reported pre- 
viously [6-8]. These demonstrated good detector effi- 
ciency and particle identification capability for events of 
this kind. 

This paper reports analyses of the 3(7r+7r~), 
2(7r"'"7r~7r'^) and K~^K~2{Tr~^Tr~) final states produced 
in conjunction with a hard photon, assumed to result 
from ISR. A clear J/tp signal is observed for each of 
these hadronic final states and the corresponding J/tp 
branching fractions are measured. While BABAR data are 
available at effective cm. energies up to 10.58 GcV, the 
present analysis is restricted to energies below 4.5 GeV 
because of the increase with energy of the backgrounds 
from non-ISR multihadron production. 



II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the 
BABAR detector at the PEP-II asymmetric e+e~ storage 
ring. The total integrated luminosity used is 232 fb~^, 
which includes data collected at the T{4S) resonance 
mass (211 fb~^), and at a cm. energy 40 MeV lower 
(21 fb"^). 

The BABAR detector is described elsewhere [9|. 
Charged particles are reconstructed in the BABAR track- 
ing system, which comprises the silicon vertex tracker 
(SVT) and the drift chamber (DCH). Separation of pi- 
ous and kaons is accomplished by means of the detector of 
internally reflected cherenkov light (DIRC) and energy- 
loss measurements in the SVT and DCH. The hard ISR 
photon and photons from tt" decays are detected in the 
electromagnetic calorimeter (EMC). Muon identification 
is provided by the Instrumented flux return (IFR). 

The initial selection of candidate events requires that 
a high-energy photon in the event with E'^ ^ > 3 GeV 
be found recoiling against six good-quality charged tracks 
with zero net charge or against four good-quality charged 
tracks with zero net charge and four or more photons 
with energy higher than 0.02 GeV. Almost every candi- 
date event has extra soft photons with energy above this 
threshold, mostly due to secondary hadron interactions 
and machine background. Each charged track is required 
to originate close to the interaction region, to have trans- 
verse momentum greater than 0.1 GeV/c and to have a 
polar angle in the laboratory frame with respect to the 
collision axis in the range from 0.4 to 2.45 radians. These 
selections guarantee the quality of the charged tracks in 
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the DCH. The charged track vertex is used as the point 
of origin to calculate the angles for all detected photons. 
Events with electrons and positrons are removed on the 
basis of associated EMC energy deposition and energy- 
loss (dE/dx) information from the DCH. 

In order to study the detector acceptance and efR- 
ciency, we developed a set of simulation programs for 
radiative processes. The simulation of the 3(7r+7r~)7 
and 2(7r+7r^7r")7 final states is based on the genera- 
tor developed according to the approach suggested by 
Kiihn and Czyz [10]. For the acceptance study we 
simulate six-charged pions in a phase space model and 
in a model which assumes only one p{770) per event, 
i.e. a /9(770)2(7r+7r-) final state. The 2(7r+7r-7r°) and 
K~^K~2{TT^Tr~) final states are simulated according to 
phase space. 

Multiple soft-photon emission from the initial-state 
charged particles is implemented with the structure- 
function technique [11, 12], while extra photon radiation 
from the final-state particles is simulated by means of 
the PHOTOS package [13]. The accuracy of the radia- 
tive corrections is about 1%. 

A sample of about 400,000 events were generated with 
these tools for each mode and passed through the de- 
tector response simulation [14]. These events were then 
reconstructed through the same software chain as the ex- 
perimental data. Variations in detector and background 
conditions were taken into account. 

For purposes of background estimation, a large sam- 
ple of events from the main ISR processes {2ttj, Sn-f ... 
57r7, 2if7r7 ...) was simulated. This sample exceeded 
the expected number of events in the dataset by a fac- 
tor of about three. In addition, the expected numbers 
of e"'"e~ <M [q — u,d,s,c) events were generated via 



JETSET [15] and e+e" 



via KORALB [16] in 



order to estimate background contributions from non- 

ISR events. The cross sections for the above processes 
are known with about 10% accuracy or better, which is 
sufficient for the background contribution study. 



III. THE KINEMATIC FIT PROCEDURE 

The initial sample of candidate events is subjected to 
a constrained kinematic fit in conjunction with charged- 
particle identification to extract events corresponding to 
the final states of interest. 

For each particular six-chargcd-particle candidate, and 
for each possible combination of particle types (i.e. 
3(7r+7r~) or K'^ K~2{TT~^Tr~)), a one-constraint (IC) kine- 
matic fit is performed without using information from 
the detected photon candidate. The only constraint used 
is zero photon mass. Due to the excellent resolution of 
the DCH, the three-momentum vector of the photon is 
better determined through momentum conservation than 
through measurement in the EMC. As a consequence, the 
calibration accuracy of the EMC and its alignment with 
respect to the DCH do not contribute to the systematic 




70, 80 

xm 

FIG. 1: The one-constraint distributions for data (un- 
shaded histogram) and MC 3(7r+7r~)7 simulation (shaded his- 
togram) for six-charged-track events fitted to the six-pion hy- 
pothesis. The cross-hatched histogram is the estimated back- 
ground contribution from non-ISR events obtained from JET- 
SET. The signal and control regions are indicated. 



uncertainties. The initial e+e^ and final-state charged- 
particle four-momenta and their covariance matrices are 
taken into account. 

The fit for the six-pion final-state hypothesis is re- 
tained for every event. If only one track is identified as 
a kaon, or if two oppositely-charged kaons are identified, 
the Ar+/f ^2(7r+7r^) fit is also retained. 

For the 2(7r+7r^7r'^)7 events a kinematic fit is per- 
formed using the initial e+e", final-state charged-particle 
and photon four-momenta and their covariance matrices. 
The highest cm. energy photon is assumed to be from 
ISR. Only the direction of the photon momentum vector 
is used in the fit, not the measured energy. All other 
photons with energies above 20 MeV are paired. Com- 
binations lying within ±35 MeV/c^ of the tt" mass are 
tested, and the event combination with the best value 
is retained, subject to the additional constraint that the 
two, two-photon pairs are consistent with the tt^ mass. 
In total five constraints (5C fit) are applied. The three- 
momentum vectors obtained from the fit for each charged 
track and photon are used in further calculations. 



IV. THE 3(7r+7r-) FINAL STATE 

A. Additional selection criteria 

The results of the IC fit to the six charged- track can- 
didates are used to make the final selection of the six- 
pion sample. The momentum vector of the photon re- 
constructed by the fit in the laboratory frame is required 
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FIG. 2: The six-pion invariant mass distribution (unshaded 
histogram) for the signal region of Fig. 1. The points indicate 
the background estimated from the difference between data 
and MC events for the control region of Fig. 1, normalized 
to the difference between data and MC events in the signal 
region of Fig. 1. The cross-hatched histogram corresponds to 
the non-ISR background of Fig. 1. 



to have a polar angle 0^^ in the range from 0.35 to 2.4 ra- 
dians and to match the measured polar angle 0™'^^ of the 
ISR photon in the EMC within 50 mrad. The correspond- 
ing azimuthal angles, (/>^* and 0™®'^'^, are also required to 
agree within this same tolerance. These angular criteria 
reduce the background by a factor of about two with no 
noticeable loss of signal. Finally, the polar angle 6*^^ of 
each charged track obtained from the fit has to satisfy 
0.45 < < 2.4 radians in order to fall within the ac- 
ceptance of the DIRC, which provides about 80% of the 
kaon identification efficiency. 

The IC-fit distribution for the six-pion candidates is 
shown as the upper histogram of Fig. 1, while the shaded 
region is for the corresponding MC-simulated pure Gn'y 
events. The experimental distribution has a contribu- 
tion from background processes, but the pure Gttj MC- 
simulated distribution is also much broader than the 
usual one-constraint distribution. This is due to mul- 
tiple soft-photon emission (detected or not detected) in 
the initial state and radiation from the final-state charged 
particles, neither of which is taken into account by the 
constrained fit but which exist both in the data and the 
MC simulation. The MC simulated distribution of 
Fig. 1 is normalized to the data in the region < 1 
where the background contamination and multiple soft- 
photon emission due to ISR or FSR is lowest. 

The cross-hatched histogram in Fig. 1 represents the 
non-ISR background contribution obtained from the 
JETSET simulation of quark-anti-quark production and 
hadronization and does not exceed 8%. 



We require xItt < 20 for the six-pion hypothesis, and 
that any accompanying fit to the 2K4:t: hypothesis have 
xIkatt > 20. The subscripts "67r" and "2is:47r" here and 
below refer to the 3(7r"'"7r~) and K~^K~2{n^TT~) final 
states exclusively. We estimate that these requirements 
reduce the misidentification of 2isr47r events from 11% to 
about 2%, at the cost of the loss of about 5% of the signal 
Gtt events. 

The region 20 < xItt < 40 is chosen as a control region 
for the estimation of background from other ISR and non- 
ISR multi-hadron reactions. The procedure followed is 
described in the next section. 

The signal region of Fig. 1 contains 19,683 data and 
19,980 MC events, while for the control region the corre- 
sponding numbers are 2,021 and 875 respectively. 



B. Background estimation 

The non-ISR background contribution to the signal re- 
gion is obtained from the JETSET MC simulation, nor- 
malized using the integrated e+e^ luminosity. The x^ 
distribution for non-ISR events is shown by the cross- 
hatched histogram of Fig. 1. The non-ISR background 
dominates by e+e~ — * dhadrons + 7r° production with 
photon from tt" mistakenly taken as ISR photon. 

MC simulation of the r+r~ final state and ISR pro- 
duction of multi-hadron final states other than 3(7r+7r~) 
shows that such states yield a background in the selected 
six-pion sample that exhibits a relatively flat contribu- 
tion to the xItt distribution. To validate these estimates 
of backgrounds with the data, we subtract the MC simu- 
lated signal distribution (the shaded histogram of Fig. 1) 
from the unshaded one, after the non-ISR background is 
subtracted. The shape of the resulting histogram is well 
described by MC simulation of remaining background 
processes. Its absolute normalisation is used to estimate 
the level of those backgrounds in the signal region. 

The background contribution to any distribution other 
than x^ is estimated as the difference between the distri- 
butions in the relevant quantity for data and MC events 
from the control region of Fig. 1 , normalized to the dif- 
ference between the number of data and MC events in 
the signal region. The non-ISR background is subtracted 
separately. 

For example, Fig. 2 shows the six-pion invariant mass 
distribution up to 4.5 GeV/c^ for the signal region of 
Fig. 1. The points with error bars show the ISR back- 
ground contribution obtained in the manner described 
from the control region of Fig. 1. The cross-hatched his- 
togram in Fig. 2 represents the non-ISR background con- 
tribution obtained from the JETSET MC simulation. 

Both backgrounds are relatively small at low mass 
(about 6-8%), but the non-ISR background accounts for 
about 20-25% of the observed data at approximately 
4GeV/c2. 

Accounting for uncertainties in cross sections for 
background processes, uncertainties in normalization of 
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events in the control region and statistical fluctuations 
in the number of simulated events, we estimate that this 
procedure for background subtraction results in a system- 
atic uncertainty of less than 3% in the number of signal 
events in the 1.6-3 GeV/c^ region of six-pion mass, but 
that it increases to 3-5% in the region above 3 GeV/c^. 

By selecting a "background-free" 6n^ sample with only 
six charged tracks and only one photon (about 5% of 
events) we can compare distributions for data and MC 
events up to =1000. We estimate that for a Xen < 
20 selection the net signal size should be increased by 
(3 ± 2)% to allow for a slight shape difference between 
the MC and experimental distributions. 
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FIG. 3: (a) The six-pion mass distributions from MC simu- 
lation for the signal (unshaded) and control (shaded) regions 
of Fig. 1. (b) The mass dependence of the net reconstruction 
and selection efficiency obtained from simulation. The curve 
is a fit described in the text. 



C. Tracking efficiency 

The procedure to measure the track-finding efficiency 
is described in our previous paper [8] for the four-pion 
final state. The method uses events that have three 
charged-particle tracks and a hard photon. These events 
are subjected to a one-constraint fit, which uses all mea- 
sured parameters of the three tracks and the photon and 
yields the three-momentum vector of the missing charged 
pion in the laboratory frame assuming this is the only 
undetected track. If the of the fit is less than 30 and 
this vector lies within the acceptance of the DCH, the 
event is included in the data sample. The ratio of thrce- 
to four-charged track events gives the track-finding effi- 
ciency. The same procedure is applied to MC-simulated 
events. The track-finding efficiency is better for MC- 
simulated events by (0.8 ± 0.5)% per track independent 



from momentum [8]. Assuming no increase in correla- 
tions due to higher multiplicity, we apply an overall cor- 
rection of + (5±3)% to the observed six-pion event sample 
based on the previous study. 



D. Detection efficiency from simulation 
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FIG. 4: (a) The distribution in track-pair opening angle for 
the minimum of the 15 values possible for each event; (b) the 
distribution in polar angle, and (c) the transverse momen- 
tum distribution for all pions from all events. All quantities 
are in the laboratory frame; the points are for data and the 
histograms are obtained from MC simulation. 

The selection procedures applied to the data are also 
applied to the MC-simulated event sample. The re- 
sulting six-pion invariant-mass distribution is shown in 
Fig. 3(a) for the signal and control (shaded histogram) 
regions. The mass dependence of the detection efficiency 
is obtained by dividing the number of reconstructed MC 
events in each 25 MeV/c^ mass interval by the number 
generated in this same interval. Note that the detection 
efficiency computed that way is insensitive to the actual 
shape of the mass distribution of Fig. 3(a) used in MC 
simulation. The result is shown in Fig. 3(b); the curve 
is obtained from a 3rd order polynomial fit to the distri- 
bution. The efficiency falls off gradually with increasing 
mass from about 20% at 1.6 GcV/c^ to about 14% at 
4.5 GeV/c^. This efficiency estimate takes into account 
the geometrical acceptance of the detector for the final- 
state photon and the charged pions, the inefficiency of 
the several detector subsystems, and event-loss due to 
additional soft-photon emission from the initial and final 
states. 

As mentioned in Sec. II, the model used in the MC 
simulation assumes that the six-pion final state results 
predominantly from the p(770)2(7r"'"7r~) production pro- 
cess. In general, this model describes well the distribu- 
tions in many of the kinematic variables characterizing 
the six-pion final state. Some examples are shown in 
Figs. 4 and 5, in which the points with error bars rep- 
resent data while the histograms are obtained from MC 
simulation. Figure 4(a) shows the distribution in ^min, 
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FIG. 6: The e^e~ cm. energy dependence of the 3(7r'''7r~) 
cross section measured with ISR data at BABAR compared 
with the direct e'^e~ measurements by DM2 and MARK-II. 
Only statistical errors are shown. 



FIG. 5: The angular distribution of the lowest-momentum 
pion (left) and of the sum of the remaining five most energetic 
pions (right) in the six-pion rest frame with respect to the 
direction of the six-pion system in the laboratory frame for 
the five regions of six-pion mass indicated in the right hand 
plots. The fourth slice is chosen to correspond to the J/tp 
region. The points are data, and the histograms are MC 
simulation. 



the minimum chargcd-pion-pair opening angle for each 
event, while Fig. 4(b) and Fig. 4(c) represent the dis- 
tribution in polar angle, 9ch, and transverse momentum, 
Pt-, respectively, for all final-state pions. All quantities 
are calculated in the laboratory frame. The overall agree- 
ment between MC simulation and data is very good. Fig- 
ure 5 compares the distributions in cos 9, where is the 
angle between a charged pion in the six-pion rest frame, 
and the direction of the six-pion system in the laboratory 
frame. Data and MC are in rather good agreement. 

In the six-pion rest frame, the angular acceptance is 
rather uniform. A simulation without resonances using 
only six-pion phase space does not produce discernible 
deviations from the observed angular distributions, and 
does not change the overall acceptance more than by 3%. 
This value is taken as an estimate of the model-dependent 
systematic uncertainty in the acceptance. 



E. Cross section for e+e — ^ 3(77"*" tt ) 

Data from the reaction e^e^ fx^fi^j are used 
to convert the invariant-mass distribution for an ISR- 
produced hadronic final state to the energy dependence 



of the corresponding e+e" cross section. The invariant 
mass of the muon pair m^^^ defines an effective e"'"e~ cm. 
collision energy, E'c.m.- The differential luminosity, dC, 
associated with the interval dEc.m. centered at effective 
collision energy Ec.m. is then obtained from 



where Er 



dN^^.,{E^ 



(l+'5FSR)-^P+p-(^c.m.)-(l+<5vac) ' 



TO 



MM . 



dN„ 



(2) 

is the number of muon 



pairs in the mass interval c?to(J 



dEr, 



is the 



acceptance, corrected for muon identification and soft- 
photon emission; (1 + (5pgj^) corrects for hard photon 
emission from final-state muons; (^-c.m.) is the 

e'^e~ — !■ fJ.^ IJ-^ Born cross section at center-of-mass en- 
ergy Ec.m.', and (1 + 5vac) is the corresponding vacuum 
polarization correction [17]. The dependence of the dif- 
ferential luminosity on Ec.m. is presented in our previous 
paper [8]. 

From a detailed study of the e^e~ fi'^iJ,~j detec- 
tion and identification efficiency described in Ref. [6] and 
comparison of the observed invariant-mass spectrum with 
theoretical calculations, we estimate the systematic un- 
certainty associated with luminosity determination to be 
3%. 

The six-pion e+e^ cross section can then be calculated 
from 



a(3(7r+7r-))(^e.„J = 



dNQ.^j{Ec.m.) 



dC{Ec 



■=671 



^6ir l-C/c.m.J 



where TO; 



6tt 



(3) 

£-0.111. with TO?^ the invariant mass of the 
six-charged-pion system; dNQ^^y is the number of selected 
six-charged-pion events after background subtraction in 
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the interval dEc.m. and e^^(iJc.m.) is the corresponding 
detection efficiency obtained from the MC simulation. 
The factor eg™'' takes into account the difference between 
the distributions for data and MC events, and the 
tracking-efficiency discrepancies discussed in Sec. IV B 
and Sec. IV C respectively. 

The energy dependence of the cross section for the re- 
action e~^e~ —^ 3(7r"'"7r~) after all corrections is shown in 
Fig. 6. It shows a structure around 1.9 GeV, reaches a 
peak value of about 1.5 nb near 2.0 GeV, followed by 
a monotonic decrease toward higher energies perturbed 
only by a peak at the J/ip mass position. The cross sec- 
tion for each 25 MeV interval is presented in Table I. 

Since dC has been corrected for vacuum polarization 
and final-state soft-photon emission, the six-pion cross 
section measured in this way includes effects due to vac- 
uum polarization and final-state soft-photon emission. 
For g — 2 calculations, vacuum polarization contributions 
should be excluded from this data. 

We studied the resolution in six-pion mass with MC 
simulation, finding that Gaussian fits of line shapes give 
mass resolutions ares that vary between 6.2 MeV/c^ in 
the 1.5-2.5 GeV/c^ mass region and 8.7 MeV/c^ in the 
2.5-3.5 GeV/c^ mass region. The observed line shape 
is not purely Gaussian mainly due to soft-photon radi- 
ation. Since the resolution varies slowly with mass, the 
data are binned in mass intervals of 25 MeV/c^ and the 
cross section has no sharp peaks (except for the J ftp re- 
gion discussed below), unfolding the resolution has little 
effect on the measured energy dependence and has not 
been performed in this analysis. For the sake of any fu- 
ture comparisons of our data with models (e.g. for the 
structure near 1.9 GeV/c^) we provide the relation be- 
tween the observed number iV^ of events in bin i and the 
"true" number of events in nearby bins: 

= el.iV^2+e2-7V^i-l-e3-iVf-|-e4.Ar,o,i+e5-7VlV2 : (4) 



G. Physics results 

The six-charged-pion cross section measured by BABAR 
can be compared with existing e+e~ measurements per- 
formed by the DM2 [18] and MARK-II [19] detectors 
(see Fig. 6). The agreement is relatively good, but the 
BABAR measurement is much more precise. The struc- 
ture around 1.9 GeV reported by both the DM2 [18] and 
FOCUS [24] experiments is clearly seen (see discussion 
later). 




m(3(7i 71 )) (GeV/c ) m(3(7i 7i )) (GeV/c ) 



FIG. 7: Invariant mass distributions for all six pion versus 
neutral four pion combinations (top) and neutral two pion 
combinations (bottom) for data (a,c) and MC simulation 
(b,d). The J/^ signal seen in the data is not included in 
the simulation. 



where is the number of events in bin i be- 

fore accounting for resolution. The coefficients 
el, ..65=0.005,0.031,0.896,0.062,0.006 are the true event 
fractions in five energy bins of 25 MeV width for the 
six-pion mass region 1.5-2.5 GeV/c^. 



F. Summary of systematic studies 

The measured six-charged-pion cross section values 
shown in Fig. 6 and summarized in Table I include only 
statistical errors. The systematic errors discussed in pre- 
vious sections are summarized in Table II, along with the 
corrections that were applied to the measurements. 

The two systematic corrections applied to the mea- 
sured cross sections sum up to -1-8% with 6-8% taken as 
a systematic uncertainty. 



Different mass combinations were studied in data and 
MC events to search for any structures or states not in- 
cluded in the simulation. Figure 7 shows the scatter- 
plots of 47r- and 27r-mass versus 67r-mass for data and 
MC events. Good agreement is seen except for narrow 
regions around the J/ip and ip{2S) masses that are not 
included in the simulation. 

In order to make a more detailed study, five intervals 
of 67r-mass are selected: (1) 1.0-1.6 GeV/c^; (2) 1.6- 
2.0 GeV/c^ (an interval with significant structure in the 
cross section, see Fig. 6); (3) 2.0-3.05 GeV/c^; (4) 3.05- 
3.15 GeV/c^ (the Jf-ip resonance region); and (5) 3.0- 
4.5 GeV/c^. Figure 8 shows the two and three pion mass 
projections, as well as the mass projections of the remain- 
ing four pions in events having a two pion combination 
within ±75 MeV/c^ of the p mass, for the five regions of 
six-pion mass just described, for both data and simula- 
tion. Background in the data is subtracted using the 
distributions as described above. 
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TABLE I: Summary of the e^e —f 3(7r"'"7r ) cross section measurement. Errors are statistical only. 



Sc.m. (GeV) 


a 


(nb) 


Sen,. (GeV) 




(nb) 


Ec.m. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


1.3125 


0.01 


± 


0.01 


2.1125 


1.36 


± 


0.12 


2.9125 


0.55 


± 


0.08 


3.7125 


0.26 


± 0.06 


1.3375 


0.01 


± 


0.01 


2.1375 


1.35 


± 


0.11 


2.9375 


0.51 


± 


0.08 


3.7375 


0.16 


± 0.05 


1.3625 


0.01 


± 


0.01 


2.1625 


1.45 


± 


0.12 


2.9625 


0.60 


± 


0.08 


3.7625 


0.18 


± 0.05 


1.3875 


0.01 


± 


0.01 


2.1875 


1.17 


± 


0.11 


2.9875 


0.68 


± 


0.08 


3.7875 


0.24 


± 0.05 


1.4125 


0.00 


± 


0.02 


2.2125 


1.43 


± 


0.12 


3.0125 


0.70 


± 


0.07 


3.8125 


0.29 


± 0.05 


1.4375 


0.02 


± 


0.01 


2.2375 


1.38 


± 


0.11 


3.0375 


0.45 


± 


0.07 


3.8375 


0.13 


± 0.05 


1.4625 


0.01 


± 


0.02 


2.2625 


1.36 


± 


0.11 


3.0625 


0.54 


± 


0.07 


3.8625 


0.21 


± 0.05 


1.4875 


0.03 


± 


0.02 


2.2875 


1.44 


± 


0.12 


3.0875 


1.87 


± 


0.10 


3.8875 


0.17 


± 0.05 


1.5125 


0.05 


± 


0.03 


2.3125 


1.40 


± 


0.11 


3.1125 


1.58 


± 


0.10 


3.9125 


0.15 


± 0.04 


1.5375 


0.10 


± 


0.03 


2.3375 


1.28 


± 


0.11 


3.1375 


0.52 


± 


0.07 


3.9375 


0.22 


± 0.05 


1.5625 


0.12 


± 


0.03 


2.3625 


1.28 


± 


0.10 


3.1625 


0.51 


± 


0.07 


3.9625 


0.20 


± 0.05 


1.5875 


0.17 


± 


0.05 


2.3875 


1.21 


± 


0.10 


3.1875 


0.55 


± 


0.06 


3.9875 


0.19 


± 0.04 


1.6125 


0.19 


± 


0.05 


2.4125 


1.38 


± 


0.11 


3.2125 


0.51 


± 


0.07 


4.0125 


0.14 


± 0.04 


1.6375 


0.24 


± 


0.06 


2.4375 


1.10 


± 


0.10 


3.2375 


0.55 


± 


0.07 


4.0375 


0.17 


± 0.04 


1.6625 


0.35 


± 


0.06 


2.4625 


1.10 


± 


0.10 


3.2625 


0.38 


± 


0.06 


4.0625 


0.17 


± 0.04 


1.6875 


0.62 


± 


0.07 


2.4875 


1.08 


± 


0.10 


3.2875 


0.40 


± 


0.06 


4.0875 


0.14 


± 0.04 


1.7125 


0.72 


± 


0.09 


2.5125 


0.92 


± 


0.10 


3.3125 


0.53 


± 


0.06 


4.1125 


0.19 


± 0.04 


1.7375 


0.98 


± 


0.09 


2.5375 


1.08 


± 


0.09 


3.3375 


0.33 


± 


0.06 


4.1375 


0.18 


± 0.04 


1.7625 


0.96 


± 


0.11 


2.5625 


1.13 


± 


0.10 


3.3625 


0.30 


± 


0.06 


4.1625 


0.13 


± 0.04 


1.7875 


1.31 


± 


0.11 


2.5875 


1.12 


± 


0.10 


3.3875 


0.37 


± 


0.06 


4.1875 


0.15 


± 0.04 


1.8125 


1.33 


± 


0.11 


2.6125 


1.10 


± 


0.10 


3.4125 


0.35 


± 


0.06 


4.2125 


0.14 


± 0.04 


1.8375 


1.44 


± 


0.12 


2.6375 


0.93 


± 


0.10 


3.4375 


0.31 


± 


0.06 


4.2375 


0.08 


± 0.04 


1.8625 


1.35 


± 


0.12 


2.6625 


1.12 


± 


0.09 


3.4625 


0.34 


± 


0.06 


4.2625 


0.13 


± 0.04 


1.8875 


1.09 


± 


0.11 


2.6875 


0.87 


± 


0.09 


3.4875 


0.39 


± 


0.05 


4.2875 


0.13 


± 0.04 


1.9125 


1.08 


± 


0.10 


2.7125 


0.94 


± 


0.09 


3.5125 


0.32 


± 


0.05 


4.3125 


0.16 


± 0.04 


1.9375 


0.91 


± 


0.10 


2.7375 


0.86 


± 


0.10 


3.5375 


0.28 


± 


0.05 


4.3375 


0.17 


± 0.04 


1.9625 


1.14 


± 


0.10 


2.7625 


0.75 


± 


0.09 


3.5625 


0.28 


± 


0.06 


4.3625 


0.08 


± 0.04 


1.9875 


1.01 


± 


0.10 


2.7875 


0.89 


± 


0.09 


3.5875 


0.22 


± 


0.06 


4.3875 


0.04 


± 0.04 


2.0125 


1.19 


± 


0.11 


2.8125 


0.91 


± 


0.09 


3.6125 


0.30 


± 


0.05 


4.4125 


0.10 


± 0.04 


2.0375 


1.54 


± 


0.11 


2.8375 


0.75 


± 


0.08 


3.6375 


0.35 


± 


0.05 


4.4375 


0.16 


± 0.04 


2.0625 


1.49 


± 


0.11 


2.8625 


0.91 


± 


0.08 


3.6625 


0.25 


± 


0.05 


4.4625 


0.11 


± 0.03 


2.0875 


1.48 


± 


0.11 


2.8875 


0.71 


± 


0.09 


3.6875 


0.41 


± 


0.05 


4.4875 


0.06 


± 0.04 



TABLE II: Summary of systematic errors for the e'^e — > 3(7r"'"7r ) cross section measurement. 



Source Correction applied Systematic error 



Luminosity from nfi'y 




3% 


MC-data difference in < 20 signal region 


+3% 


2% 


Background subtraction 




3% for men < 3.0 GeV/c^ 
5% for men > 3.0 GeV/c^ 






MC-data difference in tracking efficiency 


+5% 


3% 


Radiative corrections accuracy 




1% 


Acceptance from MC (model-dependent) 




3% 



Total (assuming addition in quadrature and no correlations) +8% 6% for men < 3.0 GeV/c^ 

8% for men > 3.0 GeV/c^ 



A simple model with only one p(770) per event is in 
excellent agreement with experimental data. No other 
significant structures are observed. A full partial wave 
analysis (PWA) would be required in order to arrive at 
a more precise interpretation of the data but the final 
state with only one p(770) per event dominates. This 
PWA requires a simultaneous analysis of the 2(7r+7r~7r°) 
final state, which is described in the next sections. 



V. THE 2(77+77-77") FINAL STATE 
A. Additional selection criteria 

The results of the five-constraint fit (see Sec. Ill) to the 

four-charged-track and five-photon candidates are used 
to make the final selection of the 2(77+77-77°) sample. The 
polar angle 6^^ of each charged track obtained from the 
fit has to satisfy 0.45 < 0^^ < 2.4 radians. To reduce 
the background from events where an energetic photon 
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FIG. 8: Invariant mass distributions for dilTerent sub-combinations of pions arranged by rows for the five different regions of 
3(7r'''7r~) mass indicated in the central column of histograms. The points (histograms) display data (simulation). The three 
pion (left column) and two pion (middle column) plots sum over all possible combinations, while the 2(7r+7r~) combinations 
(right column) include only the remaining combinations from events selected to have one 7r"'"7r~ in the p(770) region. 



is coupled with a soft photon from background to form a 
7r°, an additional selection is applied. The cosine of the 
helicity angle of the photons of the most energetic 7r°, 
measured in the tt" rest frame, is required to be less than 
0.85. We further require all photons from 7r° to have an 
energy higher than 50 MeV. 

After these additional selections we require X47r2TrO < 
40 for the six-pion hypothesis and require that none of 
the charged tracks be identified as a kaon to reduce the 
contribution from if+i^^Tr+vr^Tr^Tr" and KKsStt final 
states. 

The five-constraint-fit distribution for the four-pion 
and five-photon candidates is shown as the unshaded his- 
togram of Fig. 9, while the shaded region is for the cor- 
responding MC-simulated pure 2(7r+7r~7r°)7 events. The 
MC-simulated distribution is normalized to the data 
in the region < 10 where contamination of the back- 
ground events and multiple soft ISR and FSR is lowest. 

The cross-hatched histogram in Fig. 9 represents the 
non-ISR background contribution obtained from the 
JETSET simulation of quark-antiquark production and 
hadronization and does not exceed 10%. 

The region 40 < x\t^2'k° ^ chosen as a control 

region for the estimation of background from other ISR 



and non-ISR multihadron reactions. 

The signal region of Fig. 9 contains 35,499 data and 
6,833 MC events, while for the control region the corre- 
sponding numbers are 8,421 and 672 respectively. 



B. Background estimation 

The background subtraction procedure for the 
2(7r+7r~7r'') final state is identical to that already de- 
scribed in Sec. IV B for the 3(7r+7r^) final state, using the 
X4tt2tt° distributions shown in Fig. 9. The unshaded his- 
togram of Fig. 10 shows the 2(7r+7r^7r'^) invariant mass 
distribution for the signal region of Fig. 9. The points 
with error bars show the ISR background contribution 
obtained in the manner described from the control region 
of Fig. 9. The cross-hatched histogram in Fig. 10 rep- 
resents the non-ISR background contribution obtained 
from the JETSET MC simulation. Both backgrounds 
are relatively small at low mass (about 10%), but the 
non-ISR background accounts for about 20-25% of the 
observed data at approximately 4 GeV/c^ . 

Accounting for uncertainties in cross sections for back- 
ground processes, uncertainties in normalization of the 
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FIG. 9: The five-constraint x distributions for data (un- 
shaded histogram) and MC 2(7r^7r~7r'')7 simulation (shaded) 
for four-charged-track and five-photon events fitted to the 
six-pion hypothesis. The cross-hatched histogram is the esti- 
mated background contribution from non-ISR events obtained 
from JETSET. The signal and control regions are indicated. 



control sample, and statistical fluctuations in the num- 
ber of simulated events, we estimate that this procedure 
for background subtraction results in a systematic un- 
certainty of less than 3% in the number of signal events 
in the 1.6-3 GeV/c^ region of six-pion mass, but that 
the uncertainty increases to 5-10% in the region above 
3 GeV/c2. 

By selecting a "background-free" 2(7r+7r^7r°)7 sam- 
ple with only four-charged tracks and only five photons 
(about 5% of events) wc find that distributions for 
data and MC have similar shapes. The ratio of MC 
events selected by 'xX^^iT^n < 1000 and X47r2TrO ^ '-^^^ 
is 1.14 reflecting soft photon radiation processes. An es- 
timate for this ratio can also be made directly from the 
data by measuring the relative number of J/ip events 
over continuum for the regions noted, yielding the 
ratio 1.08 ± 0.04, in good agreement with the estimate 
from simulation. The 6% difference between the two es- 
timates is taken as the estimate of systematic error for 
the xL27rO < 40 selection. 



C. Pion- finding efficiency 

The charged-pion tracking inefficiency is corrected by 
applying a -I- (3 ± 2)% correction to the number of ob- 
served 2(7r+7r~7r'') ISR events, following the prescription 
discussed above in Sec. IV C. 

The difference in the 7r°-finding efficiencies between 
data and MC events has been studied previously using 
the 7r+7r~7r*'7 reaction [7]. To extend this study to the 




FIG. 10; The six-pion (four charged and two neutral) invari- 
ant mass distribution (unshaded histogram) for the signal re- 
gion of Fig. 9. The points indicate the background estimated 
from the difference between data and MC events for the con- 
trol region of Fig. 9, normalized to the difference between 
data and MC events in the signal region of Fig. 9. The cross- 
hatched histogram corresponds to the non-ISR background of 
Fig. 9. 

case where there are two 7r"s in the final state, a high 
statistics sample of ISR-produced W7r° 7r+7r~7r''7r'' 
events is selected using a IC fit that ignores the tt*^ from 
the Lu decay. The Tr'^-finding efficiency is then computed 
by comparing the number of events in the to peak where 
the tt" is found to the number of events where it is not. 
By comparing data and MC results, it is found that the 
7r° efficiency for simulation is 2.8±1.0±1.0% higher where 
the systematic error comes mostly from the background 
subtraction procedure. Assuming no correlation in tt*^- 
finding efficiency for two tt'^s in the event we apply -1-5.6% 
overall correction to which we assign a systematic error 
of 3%. 



D. Detection efficiency from simulation 

The detection efficiency is determined in the same 
manner described in section IV D. The simulated 
2(7r+7r^7r'^) invariant-mass distributions after selection 
are shown in Fig. 11(a) for the signal and control (shaded 
histogram) regions. The mass dependence of the detec- 
tion efficiency is obtained by dividing the number of re- 
constructed MC events in each 25 MeV/c^ mass inter- 
val by the number generated in this same interval. The 
result is shown in Fig. 11(b) and demonstrates practi- 
cally uniform efficiency versus mass. This efficiency esti- 
mate takes into account the geometrical acceptance of the 
detector for the final-state photon and the charged and 
neutral pions, the inefficiency of several detector subsys- 
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FIG. 11: (a) The six-pion (four charged and two neutral) 
mass distribution from MC simulation for the signal and con- 
trol (shaded) regions of Fig. 9. (b) The mass dependence of 
the net reconstruction and selection efficiency obtained from 
simulation. 
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Comparison of different kinematic parameters for 
events for data (points) and MC simulation (his- 



FIG 

2(77"* TT TT 

tograms): (a) momentum distributions for the fast (shaded) 
and slow (unshaded) vr^s from MC compared with the data; 
(b) combined polar angle distribution for the tt'^s; (c) com- 
bined momentum distribution for the two charged pious; (d) 
combined polar angle distribution for the two charged pious. 



terns, and event loss due to additional soft-photon emis- 
sion from the initial and final states. 

As mentioned in Sec. II, the model used in the MC sim- 
ulation for 2(7r+7r~7r°) is pure phase space. In general, 
this model describes the distributions of all the kinematic 
variables characterizing the final state well, as demon- 



strated in Fig. 12. The uniform angular acceptance in 
the rest frame of the final state pions makes the detection 
efficiency quite insensitive to the presence of intermediate 
resonance structures. This feature has also been demon- 
strated in the six charged pion analysis discussed above 
and in our earlier study of four charged pions [8] , where 
MC simulations either with intermediate resonances or 
with phase space differ by no more than 3% in detection 
efficiency. We include this 3% here as a conservative es- 
timate of the systematic error for the model dependence 
of the 2(7r"'"7r~7r°) detection efficiency. 



DM2 
t - MUPI 

i — Gamma — Gamma 
□ - MEA 
* - BoBar 




FIG. 13: The e^e~ cm. energy dependence of the 2(7r+7r~7r'^) 
cross section measured with ISR data at BABMi compared 
with the direct e^e~ measurements by detectors at ADONE 
and DCI. Only statistical errors are shown. 



E. Cross section for e~'~e — > 2(7r~'~7r tt") 



The e^e energy dependent cross section for the 
2(7r+7r^7r") final state can then be calculated from 



cr(47r27r°)(£;c 



where 



dC{Ec.m.) ■ ^4^2TrO ' ^47r27rO (-^c.m. ) 

(5) 

Ec.m. with mf^^^'^ the invariant mass of 
the six-pion system; dN^^2T!°-t is the number of selected 
six-pion events after background subtraction in the in- 
terval dEc.m. and e^jirO (-^c.m.) is the corresponding de- 
tection efficiency obtained from the MC simulation. The 



factor e';"" takes into account the differences between 
data and MC in tracking and 7r° efficiencies, as discussed 
in Sec. VC. 

The energy dependence of the cross section for the re- 
action e~^e~ 2(7r"'"7r~7r°) after all corrections is shown 
in Fig. 13. It again shows a structure around 1.9 GeV, 
reaching a peak value of about 6 nb near 2.0 GeV. The 



16 



TABLE III: Summary of the e^e — > 2(7r"'"7r tt") cross section measurement. Errors are statistical only. 



Sc.m. (GeV) 


a 


(nb) 


E,.^. (GeV) 


a 


(nb) 


E,.^. (GeV) 


a 


(nb) 


Bern. (GeV) 


a 


(nb) 


1.3125 


0.00 


± 


0.04 


2.1125 


6.44 


± 


0.43 


2.9125 


2.46 


± 


0.26 


3.7125 


0.74 


± 0.15 


1.3375 


0.00 


± 


0.05 


2.1375 


5.45 


± 


0.42 


2.9375 


2.49 


± 


0.26 


3.7375 


0.51 


± 0.13 


1.3625 


0.04 


± 


0.06 


2.1625 


5.92 


± 


0.41 


2.9625 


2.16 


± 


0.26 


3.7625 


0.60 


± 0.14 


1.3875 


0.06 




0.07 


2.1875 


4.81 




0.40 


2.9875 


2.64 




0.26 


3.7875 


0.64 


± 0.13 


1.4125 


0.09 


± 


0.08 


2.2125 


5.81 




0.40 


3.0125 


2.88 




0.26 


3.8125 


0.86 


± 0.14 


1.4375 


0.00 


± 


0.06 


2.2375 


5.45 


± 


0.39 


3.0375 


2.01 


± 


0.26 


3.8375 


0.73 


± 0.13 


1.4625 


0.19 


± 


0.09 


2.2625 


5.68 


± 


0.39 


3.0625 


2.17 


± 


0.24 


3.8625 


0.53 


± 0.13 


1.4875 


0.04 


± 


0.09 


2.2875 


5.04 


± 


0.39 


3.0875 


6.94 


± 


0.34 


3.8875 


0.53 


± 0.12 


1.5125 


0.18 


± 


0.13 


2.3125 


5.10 


± 


0.36 


3.1125 


5.48 




0.31 


3.9125 


0.70 


± 0.13 


1.5375 


0.30 


± 


0.15 


2.3375 


4.99 


± 


0.37 


3.1375 


2.41 




0.23 


3.9375 


0.62 


± 0.13 


1.5625 


0.58 


± 


0.16 


2.3625 


5.01 


± 


0.37 


3.1625 


1.77 




0.22 


3.9625 


0.50 


± 0.12 


1.5875 


1.29 


± 


0.22 


2.3875 


4.88 


± 


0.36 


3.1875 


2.01 


± 


0.21 


3.9875 


0.55 


± 0.12 


1.6125 


2.10 


± 


0.28 


2.4125 


5.86 


± 


0.37 


3.2125 


1.63 


± 


0.20 


4.0125 


0.47 


± 0.12 


1.6375 


1.72 


± 


0.32 


2.4375 


4.07 


± 


0.33 


3.2375 


1.79 


± 


0.21 


4.0375 


0.42 


± 0.12 


1.6625 


2.88 


± 


0.35 


2.4625 


4.70 


± 


0.34 


3.2625 


1.52 


± 


0.18 


4.0625 


0.60 


± 0.12 


1.6875 


4.12 


± 


0.39 


2.4875 


4.37 


± 


0.33 


3.2875 


1.77 


± 


0.19 


4.0875 


0.40 


± 0.11 


1.7125 


4.67 




0.41 


2.5125 


4.25 




0.34 


3.3125 


1.54 


± 


0.19 


4.1125 


0.50 


± 0.11 


1.7375 


5.42 


± 


0.45 


2.5375 


4.39 


± 


0.33 


3.3375 


1.50 


± 


0.18 


4.1375 


0.46 


± 0.11 


1.7625 


4.86 


± 


0.41 


2.5625 


3.85 


± 


0.31 


3.3625 


1.40 


± 


0.18 


4.1625 


0.32 


± 0.11 


1.7875 


4.65 


± 


0.43 


2.5875 


4.13 


± 


0.31 


3.3875 


1.38 


± 


0.17 


4.1875 


0.45 


± 0.10 


1.8125 


4.97 




0.44 


2.6125 


3.92 


± 


0.31 


3.4125 


1.46 




0.17 


4.2125 


0.17 


± 0.10 


1.8375 


5.65 




0.44 


2.6375 


3.49 




0.31 


3.4375 


1.27 




0.18 


4.2375 


0.50 


± 0.11 


1.8625 


4.86 


± 


0.41 


2.6625 


3.74 


± 


0.31 


3.4625 


1.14 


± 


0.17 


4.2625 


0.14 


± 0.10 


1.8875 


4.37 


± 


0.39 


2.6875 


3.41 


± 


0.30 


3.4875 


1.35 


± 


0.17 


4.2875 


0.39 


± 0.11 


1.9125 


4.05 


± 


0.38 


2.7125 


3.72 


± 


0.30 


3.5125 


1.06 


± 


0.17 


4.3125 


0.26 


± 0.10 


1.9375 


4.44 


± 


0.39 


2.7375 


2.94 


± 


0.29 


3.5375 


1.07 


± 


0.16 


4.3375 


0.56 


± 0.10 


1.9625 


4.48 


± 


0.38 


2.7625 


3.34 


± 


0.30 


3.5625 


0.87 


± 


0.15 


4.3625 


0.26 


± 0.11 


1.9875 


4.85 


± 


0.40 


2.7875 


3.11 


± 


0.29 


3.5875 


0.70 


± 


0.15 


4.3875 


0.30 


± 0.10 


2.0125 


5.59 


± 


0.41 


2.8125 


3.43 


± 


0.29 


3.6125 


0.70 


± 


0.16 


4.4125 


0.34 


± 0.10 


2.0375 


5.15 


± 


0.40 


2.8375 


2.81 


± 


0.28 


3.6375 


1.18 


± 


0.16 


4.4375 


0.16 


± 0.10 


2.0625 


6.12 


± 


0.42 


2.8625 


3.12 


± 


0.27 


3.6625 


0.74 


± 


0.15 


4.4625 


0.32 


± 0.10 


2.0875 


5.97 


± 


0.42 


2.8875 


2.45 


± 


0.26 


3.6875 


1.65 


± 


0.17 


4.4875 


0.17 


± 0.10 



cross section for each 25 MeV interval is presented in 
Table III. 

Since dL (see eq. 2) has been corrected for vacuum po- 
larization and final-state photon emission, the cross sec- 
tion includes effects due to vacuum polarization and final- 
state soft-photon emission. For g—2 calculations, vacuum 
polarization contributions should be excluded from our 
data. 

The observed line shape is not purely Gaussian mainly 
due to soft-photon radiation. Once again, no unfolding 
of the resolution is performed for the results shown here. 
The coefficients in equation 4 appropriate to the 2(pi-|-pi- 
piO) case are el, ..e5=0.007,0.091,0.744,0.114,0.011. 



F. Summary of systematic studies 

The measured cross sections for the 2(7r+7r~7r'^) final 
state, shown in Fig. 13 and tabulated in Table III, include 
only statistical errors. The systematic errors discussed in 
previous sections are summarized in Table IV, along with 
the corrections that were applied to the measurements. 
The two systematic corrections applied to the measured 



cross sections sum to -|-8.8% with 10-13% systematic un- 
certainty. 



G. Physics results 

The cross section for the 2(7r+7r~7r°) final state mea- 
sured by BABAR (Fig. 13) can be compared with existing 
e+e" measurements performed by the /utt [20], Gamma- 
Gamma [21] and MEA [22] detectors at the ADONE col- 
lider and by the DM2 [18] detector at DCI. The BABAR 
measurement is much more precise and disagrees with 
DM2. The latter have probably large systematic errors 
due to normalization, not discussed in the original refer- 
ence [18] that mostly focuses on the confirmation of the 
dip at 1.9 GeV previously observed in the six-charged- 
pion mode. 

Different mass combinations were studied in data and 
MC events to search for any structures or states not in- 
cluded in the simulation. 

Figure 14 shows invariant mass scatter plots for; (a) 

neutral three- versus six-pion combinations, and (b) neu- 
tral 7r"'"7r~7r'^ combinations with the lower momentum tt'* 
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TABLE IV: Summary of systematic errors for the e'^e — > 2(7r+7r tv°) cross section measurement. 



Source 


Correction applied 


Systematic error 




Luminosity from fifij 


- 


3% 




MC-data difference in < 40 signal region 


0% 


6% 




Background subtraction 




5% for m2(^+^-.^o) 


< 3.0 GeV/c^ 








, > 3 CeV/r^ 


MC-data difference in tracking efficiency 


+3% 


2% 




MC-data difference in vr" losses 


+5.6% 


3% 




Radiative corrections accuracy 




1% 




Acceptance from MC (model-dependent) 




3% 




Total(assuming addition in quadrature and no correlations) 


+8.8% 


10% for m2(^+^-^o 
13% for m2(^+^-^o 


) < 3.0 GeV/c^ 
) > 3.0 GeV/c^ 




12 3. 
m(4jr2jt°) (GeV/c^) 



m(2jt7i,o) (GeV/cO 



FIG. 14: Invariant mass scatter plots from data for; (a) the 
six pion (2(7r+7r^7r")) final state versus the neutral three-pion 
combinations (7r"''7r~7r''); and (b) one neutral three-pion com- 
bination with the lower momentum tt" versus the other with 
the higher momentum tt". 



versus three-pion with the higher momentum tt*^. The 
a;(782) and 77 mesons are seen in the three body combina- 
tions. Figure 15(a) shows the projection of the three-pion 
invariant mass distribution of Fig. 14 with clear ijj(782) 
and -q signals. Figures 15(b,c,d) show mass projections 
for the two-pion combinations. Backgrounds are sub- 
tracted from the data points shown using the con- 
trol region, and the non-lSR JETSET simulation (shaded 
histograms in Fig. 15), as described above. The /9(770) 
meson is clearly seen in the tt^tt^ and tt+tt^ combina- 
tions, while there is some indication of a small signal 
from /o(980) in the tt^tt" mass distribution. The phase 
space MC simulation shown by unshaded histograms in 
Fig. 15 does not include any of these structures. 

Figure 16 compares the tt+tt^ mass distribution (mul- 
tiplied by two) with the summed tt^tt*^ mass distribu- 
tion. The basic shapes are quite similar, although there 
are more than twice as many charged as neutral p(770) 
mesons. The neutral 7r'*"7r~ also has a broad bump 
around 1.3 GeV/c^ relative to the charged tt^tt^ distri- 
bution, perhaps indicating the presence of some interme- 
diate /o(1370) or /2(1270) production. 

No obvious structures are seen in the four or five pion 
combinations (not shown). Though the structures ob- 
served above are suggestive, a partial wave analysis would 
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FIG. 15: Invariant mass distributions from data (points), 
simulation (unshaded histograms), and the non-ISR con- 
tributions obtained from JETSET (shaded histograms) for; 
(a) neutral three pion (7r^7r~7r°) combinations; (b) charged 
two pion (tt^tt") combinations; (c) neutral two body pi-zero 
(tt^tt") combinations; and (d) neutral two body charged pion 
{■k'^tx~) combinations. 



be needed to interpret the data more completely. 

Estimates for the ti;(782) and 77 contributions can 
be obtained using the ■k^i:~'t:^ mass distribution of 
Fig. 17(a) which shows an expanded view of Fig. 15(a). 
The procedure is to fit each signal with a two-Gaussian 
function representing the signal plus a polynomial back- 
ground, yielding 9899 ± 158 and 770 ±40 lott+tt-tt^ and 
TyTr+Tr^Tr*^ events, respectively. 

The energy dependence of the wtt+tt^tt'' cross section 
is determined by performing this fit for each 25 MeV/c^ 
bin of the six-pion (2(7r+7r^7r°)) mass (Fig. 18). In addi- 
tion to the broad structure at low mass, there is also 
a sharp structure in the J /i}) region corresponding to 
170 ± 24 events decaying into the a;7r+7r~7r° final state. 
After correcting for efficiency and normalizing to the ISR 
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FIG. 16: The scaled tt^tt^ mass distribution (points, with 
errors, multiplied by two) compared with the summed tt^tt" 
mass distribution (histogram). 
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FIG. 18: The cross section versus Ec.m. for all 2(7r^7r~7r'') 
events (circles), lutv^ tv~ tt'^ events (squares), and oji] events 
(triangles). 
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FIG. 17: The 7r+7r~7r° invariant mass distributions; (a) in 
the low mass region; and (b) for the second 7r"'"7r~7r'' invari- 
ant mass combination remaining in those events where the 
first combination lies within the mass range indicated by the 



arrows m a 



luminosity, this yields the e^e" — s- wTr+Tr^Tr*^ cross sec- 
tion shown by the squares in Fig. 18. 

Six pion events that contain an u, defined by the ar- 
rows in Fig. 17(a), also sometimes contain an rj as shown 
in Fig. 17(b). After selecting these uirj events, and sub- 
tracting the background using the 77 side bands, we calcu- 
late the cross section for e^e~ — > ur] presented in Fig. 19. 
The cross section is corrected for the decay rate of w and i] 
to TT+Tr^Tr" taken from PDG [23]. A prominent structure 
can be seen around 1.6 GeV, with a smaller peak (from 
13 events over a background of less than 0.5) around 3.1 
GeV from J/tp ^ ujrj decay. The observed cross section 
is fitted with a resonance-type parameterization 



^m(s) = 



■m,^/'^^o^y (Jo/ F{m'^) 



(6) 
(7) 



where m is the mass, Fq the width, and (Tq the peak cross 
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FIG. 19: The ujri cross section from the 2{tt^-k~ tt'^) event 
sample. The line is the fit to the structure in the 1.6 GeV 
region described in the text. 



section of this resonance production in e+e^ collisions. 
F{s) is a phase space term equal to the cube of the lo 
(or 77) momentum in the Lurj rest frame. The fit gives 
TO = 1.645 ±0.008 GeV/c2, Fq = 0.114 ±0.014 GeV, and 
(To = 3.08 ± 0.33 nb. 

As discussed in section VE, the effect of mass resolu- 
tion on the measured width of this structure is small and 
not taken into account. The mass value obtained is close 
to the value 1670 ± 30 GeV/c^ hsted in the PDG [23] 
for the a;(1650), for which oj-q final state has been seen 
by earlier experiments, but the width seen here is sub- 
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FIG. 20: Ratio of the 2(7r^7r~7r") to 3(7r"^7r~) cross sections 
for; (a) selected events; (b) the same events excluding the 
LUTv'^Tv~ Tv^ contribution from the 2(7r'''7r~7r°) sample. The 
lines are fits described in the text. 
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FIG. 21: The 3(7r^7r~) cross section data, as given in Fig. 6, 
compared with the coherent fit (sofid line) between resonance 
and continuum terms described in the text. The dashed line 
and the shaded function show the individual incoherent con- 
tributions from continuum and resonance terms. 



stantially narrower than the 0.315 ±0.035 GeV listed. A 
structure decaying to lut] in this region might also corre- 
spond to the (/)(1680), whose mass and width are hsted in 
the PDG as 1.680±0.020 GeV/c^, and 0.150±0.050 GeV, 
respectively. However, no branching fraction of 0(1680) 
to the LOT] final state was reported in the previous exper- 
iments. 

Figure 20(a) shows the cross section ratio for 



2(7^^ 



to 3(7r+7r ), as given by Figs. 18 and 6, 



respectively. A good fit to this ratio for all energy inter- 
vals, shown by the curve in Fig. 20(a), is obtained with 
a constant plus a Gaussian in the 1.6 GeV region. The 
ratio equals 3.98 ± 0.06 ± 0.41 everywhere, except for the 
region around 1.6 GeV, where it reaches about 8 at peak. 
The structure may be at least partially explained by the 
presence of the ojrj structure in the 2(7r+7r^7r°) final state. 
When the ujTr~^TT~Tr^ contribution is subtracted from the 
2(7r+7r~7r°) final state, the ratio is no longer fiat, as shown 
by Fig. 20(b), and the structure at low mass disappears. 
A 2nd-order polynomial fits the data well as shown by 
the curve in Fig. 20(b). 

Figures 21 and 22 compare the 3(7r+7r^) and 



2(7r^ 



cross section data, respectively, with fits to 



the model presented in Ref. [24] . The structures observed 
in both channels around 1.9 GeV are not well described 
by a single Breit-Wigner resonance, and may result from 
rather complicated physics, such as several vector states 
decaying to the same mode. The model [25] has the form 



s — + i-y/sF 



(8) 



where m, F and (j) ^re the mass, width and relative 
phase of the Breit-Wigner type amplitude, represent- 



+ 



1 1 1 1 




\ 

1 : 




I'm 

h 


yi 






\ 






v. ■ 






''^^U 



1.5 



2.5 3 3.5 4 4.5 
E_ (GeV) 



FIG. 22: The 2(7r"'"7r~7r'') cross section data, as given in 
Fig. 18, compared with the coherent fit (solid line) between 
resonance and continuum terms described in the text. The 
dashed line and the shaded function show the individual in- 
coherent contributions from continuum and resonance terms. 



ing the structure, and 17 is a coupling constant and 
Acont = Co + ci |^^_^°^-)2-a is a Jacob-Slansky am- 
plitude [26] representing an amalgamation of broad 
resonances with co,ci,a,6, toq free parameters. The 
following "resonance" parameters are obtained for the 
structure: 
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me-n = 1.88±0.03GeV/c2;m4^27rO = 1.86±0.02 GeV/c^, 
Tq„ = 0.13 ± 0.03 GeV; r4^27r0 = 0.16 ± 0.02 GeV, 

ct)6^ = 21 ± 40° ; <Pi^2^o = -3 ± 15° . 

The parameter values obtained seem to be essentially 
independent of the final state charge combination. These 
values may also be compared with those obtained with 
a similar model by the FOCUS experiment [24, 27] m = 
1.91±0.01 GeV/c2,r = 0.037± 0.013 GeV, ^ = 10±30° . 
The mass values arc consistent, but the widths obtained 
by BABAR are substantially larger. 




X(2K43i) 



FIG. 23: The one-constraint x distributions for data (un- 
shaded histogram) and MC K'^ K~2{Tr'^Tr~)j simulation 
(shaded histogram) for six-charged-track events fitted to the 
K'^ K~2{n'^Tr~) hypothesis. The cross-hatched histogram is 
the estimated background contribution from non-ISR events 
obtained from JETSET. The signal and control regions are 
indicated. 



VI. THE ii:+Ji£:-2(7r+7r-) FINAL STATE 

The constrained fit of the six-charged-track events to 

the hypothesis of two oppositely charged kaons and four 
charged pions, where at least one of the kaons has positive 
particle identification, allows us to select this final state. 
Figure 23 shows the distributions for both data and 
simulation, where the simulation of the K~^K~ 2{7r'^7r~) 
reaction uses a phase space model with a cross section 
energy dependence close to that which we observe exper- 
imentally, and ISR and FSR extra radiative processes are 
included. Also shown is the estimated contribution from 
non-ISR events obtained by JETSET simulation. 

Figure 24(a) presents the simulated mass distribution 
for the 2K4tt events. The mass dependence of the ef- 
ficiency, calculated as a ratio of selected to generated 
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FIG. 24: (a) The K~2{tt~^tt~) mass distributions from 
MC simulation for the signal (unshaded) and control (shaded) 
regions of Fig. 23. (b) The mass dependence of the net recon- 
struction and selection efficiency obtained from simulation. 
The curve is a fit described in the text. 



2KAtt MC events, is shown in Fig. 24(b). The efficiency 
falls gradually from about 15% at low mass, to about 
3% at 4.5 GcV/c^. The data are well represented by the 
3rd-order polynomial shown by the curve in Fig. 24(b). 

After the selections X^katt < ^0, Xgtt > 20 arc im- 
posed, we estimate from simulation that about 0.6% of 
the six charged pion events remain in the sample due to 
misidentification of pions. Although the background sub- 
traction procedure, as described above, using the control 
region 20 < X2Ki-jT < ^0 is invoked, background from 
pion misidentification will not be subtracted by this pro- 
cedure. Thus, 0.6% of the six pion events shown in Fig. 2 
are subtracted, leading to a correction of about 3% to the 
total number of K^2{'k^tt^) events. 

Figure 25 shows the K^2{tt^tt^) invariant mass 
distributions for signal events, selected as defined above, 
as well as for events from the control region in X2K4tt' 
and for the expected non-ISR background from JETSET 
MC simulation. These latter two distributions are used 
to subtract background in the signal region. Clear J/tjj 
and tp{2S) signals are seen. 

Using the number of observed c;vents, efficiency, and 
ISR luminosity, we obtain the e+e^ K'^K~2{TT~^n~) 
cross section shown in Fig. 26. No measurements are 
available from earlier experiments. Table V presents the 
cross section in 25 MeV bins. The systematic errors are 
dominated by the uncertainty in the acceptance simula- 
tion (10%), by the uncertainty in the background sub- 
traction (5%), and by the difference between the kaon 
identification efficiencies for data and MC events (up to 
2% per track), and are estimated to be about 15%. 

Figures 27(a,b) show the Kn mass combinations. Pro- 
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TABLE V: Summary of the e^e — » K~^K 2{iv~^tt ) cross section measurement. Errors are statistical only. 



Sc.m. (GeV) 


a 


(nb) 


Sen,. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


Ec.m. (GeV) 


a 


(nb) 


2.0125 


0.01 


± 


0.01 


2.6375 


0.29 


± 


0.08 


3.2625 


0.36 


± 


0.10 


3.8875 


0.31 


± 0.11 


2.0375 


0.01 


± 


0.01 


2.6625 


0.44 


± 


0.10 


3.2875 


0.54 


± 


0.10 


3.9125 


0.21 


± 0.10 


2.0625 


0.00 


± 


0.00 


2.6875 


0.55 


± 


0.10 


3.3125 


0.43 


± 


0.10 


3.9375 


0.33 


± 0.09 


2.0875 


0.01 


± 


0.01 


2.7125 


0.28 


± 


0.09 


3.3375 


0.14 


± 


0.10 


3.9625 


0.25 


± 0.10 


2.1125 


0.02 


± 


0.02 


2.7375 


0.47 


± 


0.10 


3.3625 


0.38 


± 


0.10 


3.9875 


0.11 


± 0.09 


2.1375 


0.02 


± 


0.03 


2.7625 


0.29 


± 


0.09 


3.3875 


0.39 


± 


0.10 


4.0125 


0.22 


± 0.09 


2.1625 


0.08 


± 


0.04 


2.7875 


0.54 


± 


0.10 


3.4125 


0.17 


± 


0.10 


4.0375 


0.20 


± 0.11 


2.1875 


0.05 


± 


0.03 


2.8125 


0.53 


± 


0.11 


3.4375 


0.43 


± 


0.11 


4.0625 


0.17 


± 0.10 


2.2125 


0.12 


± 


0.04 


2.8375 


0.49 


± 


0.10 


3.4625 


0.34 


± 


0.10 


4.0875 


0.16 


± 0.10 


2.2375 


0.10 


± 


0.03 


2.8625 


0.32 


± 


0.10 


3.4875 


0.40 


± 


0.11 


4.1125 


0.18 


± 0.10 


2.2625 


0.02 


± 


0.03 


2.8875 


0.45 


± 


0.11 


3.5125 


0.29 


± 


0.10 


4.1375 


0.40 


± 0.10 


2.2875 


0.18 


± 


0.06 


2.9125 


0.25 


± 


0.10 


3.5375 


0.10 


± 


0.11 


4.1625 


0.21 


± 0.08 


2.3125 


0.15 


± 


0.05 


2.9375 


0.68 


± 


0.11 


3.5625 


0.17 


± 


0.10 


4.1875 


0.12 


± 0.09 


2.3375 


0.21 


± 


0.06 


2.9625 


0.42 


± 


0.10 


3.5875 


0.17 


± 


0.09 


4.2125 


0.25 


± 0.09 


2.3625 


0.08 


± 


0.05 


2.9875 


0.52 


± 


0.11 


3.6125 


0.18 


± 


0.09 


4.2375 


0.11 


± 0.09 


2.3875 


0.25 


± 


0.07 


3.0125 


0.34 


± 


0.10 


3.6375 


0.06 


± 


0.09 


4.2625 


0.06 


± 0.09 


2.4125 


0.34 


± 


0.08 


3.0375 


0.43 


± 


0.11 


3.6625 


0.24 


± 


0.09 


4.2875 


0.20 


± 0.08 


2.4375 


0.39 


± 


0.08 


3.0625 


0.69 


± 


0.11 


3.6875 


0.67 


± 


0.12 


4.3125 


0.17 


± 0.09 


2.4625 


0.38 


± 


0.08 


3.0875 


2.03 


± 


0.20 


3.7125 


0.30 


± 


0.10 


4.3375 


0.34 


± 0.10 


2.4875 


0.19 


± 


0.07 


3.1125 


1.14 


± 


0.16 


3.7375 


0.30 


± 


0.10 


4.3625 


0.20 


± 0.11 


2.5125 


0.36 


± 


0.08 


3.1375 


0.42 


± 


0.13 


3.7625 


0.18 


± 


0.11 


4.3875 


0.04 


± 0.10 


2.5375 


0.24 


± 


0.08 


3.1625 


0.62 


± 


0.12 


3.7875 


0.37 


± 


0.10 


4.4125 


0.14 


± 0.10 


2.5625 


0.30 


± 


0.08 


3.1875 


0.39 


± 


0.11 


3.8125 


0.25 


± 


0.09 


4.4375 


0.24 


± 0.09 


2.5875 


0.32 


± 


0.08 


3.2125 


0.56 


± 


0.11 


3.8375 


0.19 


± 


0.09 


4.4625 


0.31 


± 0.11 


2.6125 


0.26 


± 


0.08 


3.2375 


0.37 


± 


0.10 


3.8625 


0.22 


± 


0.10 


4.4875 


0.26 


± 0.10 



duction of K-k pairs is dominated by the K*^(Wi2) clearly 
seen in Fig. 27(a) and the projection plot in Fig. 27(b). 
Some small structure is also seen at about 1.9 GeV/c^. 
Figure 27(c) shows the K^K~ mass distribution. No 
structures except a smah signal from 0(1020) are seen. 
These events correspond to J/tp (/)(1020)2(7r+7r~) de- 
cay. This is demonstrated in Fig. 27(d) where events 
from the 1.01-1.03 GeV/c^ mass interval of Fig. 27(c) 
arc shown. We find 35 ± 7 events in the J/ijj peak from 
the above decay chain. 



culate the products: 



N{J/tlj 3(7r+7r-)) 

dC/ dE ■ CMC 
57.1 ±3.8 ±3.4 nbMeV , 
iV(J/V^3(7r+^-))-m2/^ 



Q-k'^ ■ dC/dE- CMC -C 
= (2.37 ± 0.16 ± 0.14) X 10"^ keV 



VII. THE J/V REGION 

Figure 28 shows an expanded view of the J /ijj mass re- 
gion in Fig. 2 for the six-pion data sample with no back- 
ground subtraction. The signals from J/ip ^ 3(7r+7r~) 
and i}{2S) 3(7r+7r~) are clearly seen. The obser- 
vation of J/tjj decaying into six-charged pions gives a 
direct measurement of the 67r-niass resolution and the 
absolute energy scale. A fit with a Gaussian for the J/tp 
peak and a polynomial function for the continuum gives 
CmCSTr) =8.7 McV/c^ and less than 1 MeV/c^ difference 
from the PDG [23] value for the J/tp mass. The observed 
mass resolution agrees with the simulation within 10%. 

The observed 496 ± 33 events at the J/tp peak can 
be used to calculate the branching fraction for J/tp 
3(7r+7r~). The simulation shows that because of radia- 
tive effects only 90% of the signal events are under the 
Gaussian curve. Using the corrected number, we can cal- 



where 

CT;-^// = 67r2r^/'^ C/m^/^ = 12983 ± 409 nb - MeV 

is the integral over the J/tp excitation curve; is the 
electronic width; dC/dE = 65.6 nb"VMeV is the ISR 
luminosity at the J/tjj mass; cmc = 0.147 is the detection 
efficiency from simulation with the corrections and errors 
discussed in Sec. IV F; and C = 3.894 x lO" nbMeV^ is 
a conversion constant. The subscript "67r" for branching 
fractions refers to the 3(7r+7r^) final state exclusively. 

Using Fe/'^ = 5.40 ± 0.18 keV [23], we obtain 
= (4.40 ± 0.29 ± 0.29) x 10"'', substantially 
more precise than the current PDG value Bj/^^q^ = 
(4.0 ± 2.0) X 10~^ [23]. The systematic error includes a 
3% uncertainty in . 

The tp(2S) peak mostly corresponds to the decay chain 
tp{2S) J/tpn^n^ 3(7r+7r~). The number of events 
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FIG. 25: The K^2{tv'^tv~) invariant mass distribution 
(unshaded histogram) for the signal region of Fig. 23, after 
subtracting mis-identified six charged pion events as described 
in the text. The points indicate the background estimated 
from the difference between data and MC events for the con- 
trol region of Fig. 23, normalized to the difference between 
data and MC events in the signal region of Fig. 23. The cross- 
hatched histogram corresponds to the non-ISR background of 
Fig. 23. 



extracted from a fit to a Gaussian distribution for the 
ip{'^S) peak and a polynomial function for the contin- 
uum is 61 ± 16. The direct decay of V'(2S') 3(7r+7r-) 
is very small and using -Bi/,(2S)^67r = (1-5 ± 1.0) x 10"'' 
[23] only 5 events are expected - less than the expected 
statistical fluctuation in the observed number of il>{2S) 
events. On the other hand using B^,^2S)^J/tim+-n- = 
0.317±0.011 [23] and the recently measured i?,7/0^4,r = 
(3.61 ± 0.37) X 10"-^ [8], we estimate that 42 ± -5 events 
should be seen from the tp{2S) J/^tt+tt" 3(7r"'"7r~) 
decay chain. This estimate may be compared with a di- 
rect measurement of 46 ± 8 events, using data selected 
by taking events with four-pion invariant mass combina- 
tions within ±50 MeV/c^ of the J/tp mass and counting 
the number of ijj{2S) events in the six-pion mass distribu- 
tion. Since the total number of events from the six-pion 
decay of the ip{2S) is less than la larger than the number 
decaying through the J/i/jTr~^TT~ decay chain, no signifi- 
cant measurement can be made for the direct decay rate 
of V'(25') — > 3(7r+7r~) from these data. 

Figure 29 shows the J ftp and ip{2S) signals in the 
K~^K~2{'K~^n~) mode. Background is not subtracted. 
The numbers of events under the Gaussian curves are 
205 ± 17 and 51 ± 11 respectively. As shown in Sec. VI, 
35 ± 7 J ftp events decay to the (pAn final state. The 
mass resolution is about 6.5 MeV/c^ for the 2X4^ chan- 
nel. Using B^,^2S)^J/4m+TT- ^nd the recently measured 
Bj/^^2K2tz = (6.2 ± 0.7) X 10^^ [8] we can expect 
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FIG. 26: The cm. energy dependence of the e^e~ 

K~2(n^n~) cross section obtained from ISR events at 
BABAR. Only statistical errors are shown. 
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FIG. 27: Invariant mass plots for various selections of data 

from the K'^ K~2{-k'^it^) sample for: (a) the scatter-plot of 
the K^Tv^ and K^tv^ invariant mass values; (b) the K^tx~ 
or K^Tx^ mass projection of (a); (c) the mass distribution for 
K^K~; (d) the K~2{'k~^'k~) mass distribution for events 
around the ^(1020) peak from (c). 



25 ± 3 events due the decay chain '0(25') J/tpn^TT 
K+K-2{TT+Tr-). The difference of 26 ± 13 events be- 
tween the total number of ip{2S) iC+if ~2(7r+7r~) de- 
cays and this specific decay chain provides an estimate 
of direct ^(25) ^ K+K-2{Tr+TT-) decay. 

Using the radiative correction factor 0.9, ISR luminos- 
ity dC/dE = 84.0 nb^V^^eV at ^(25) mass, and e2KA7, 
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FIG. 28: The 3(7r^7r~) mass distribution for ISR-produced 
e"'"e~ — > 3(7r"'"7r~) events in the mass region around the J/tp 
and tp{2S); there are clear signals at the J/tp and ip{2S) mass 
positions. 
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FIG. 30: The 2(7r+7r"7r°) mass distribution for ISR-produced 
events in the mass region around the J/ip and %I}{2S); there 
are clear signals at the J/ip and 'ip{2S) mass positions. 
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FIG. 29: The K+ K~2{tt+-k~) mass distribution for ISR- 
produced events in the mass region around the J/ip and 
%I}{2S); there are clear signals at the J/ip and ip{2S) mass 
positions. The latter is dominated by ip{2S) — > J / •^iir'^ 'K~ , 
with J/xp K+K'tt^tt'. 




from simulation with efRciency corrections, we obtain 

^J/rP^K+K-2{n+T7-) ' O^int 

= (66.2 ±5.5 ±4.0) nbMeV , 

= (2.75 ± 0.23 ± 0.17) x lO'^ keV , 

= (0.47 ± 0.09 ± 0.04) X lO'^ keV , 

-^K+K 2(7r+7r ) * -^ee 

= (4.4 ± 2.1 ± 0.3) X 10"^ keV . 

The systematic errors are mainly due to the uncertainties 
in the acceptance and ISR luminosity. 

Using the world average values for , F^i^'^'' and 
B^^fc+K- [23], we calculate the branching fractions 



B 



j/^^2K4. = (5.09 ± 0.42 ± 0.35) X 10-^ 
j/^^^in = (1.77 ± 0.35 ± 0.12) X 10-^ 



B 



Bi,(2S)^2Kin = (2.1±1.0±0.2) X 10-^ 

to be compared with the current world average val- 
ues [23] of (3.1 ± 1.3) X 10-3 for J/^ 2X4^ and 
(1.60 ±0.32) X 10-3 for J/tp (piir. No earlier measure- 
ments are available for tlj{2S) K~^K~2{Tr~^7T~) decays. 
The uncertainty in tU"^ has been added in quadrature 
to the systematic error estimate. 

The observation of J/ip decaying into four-charged 
and two-neutral-pions gives a direct measurement of 
the 2(7r+7r-7r°)-mass resolution and the absolute en- 
ergy scale in this mode. The mass distribution af- 
ter background subtraction is shown in Fig. 30. A fit 
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with a Gaussian for the J /ip peak and a polynomial 
function for the continuum gives a mass resolution of 
(^m(2(-K+-K--K°)) = 12.7 MeV/c^. The central value of the 
mass lies 1.5 MeV/c^ above the PDG [23] value for the 
J/ ^ mass and gives an estimate of systematic uncertainty 
in absolute energy scale determined by invariant mass of 
six pions. The observed mass resolution agrees with the 
simulation very well. 

The 761 ± 42 events observed at the ,//(,/; peak can 
be used to calculate the branching fraction for J/^ — > 
2(7r"'"7r~7r'^) as was done for J/V' — > 3(7r+7r^). We calcu- 
late the products: 



J/-4, 



N{J/i} 47r27r°) 

dC/dE ■ fMC 
215 ±12 ±24 nbMeV 
N{J/tfj 47r27rO) • m 



67r2 • dC/dE -euc-C 
(8.9 ± 0.5 ± 1.0) X 10"^ keV 



Using dC/dE = 84.0 nb"VMeV as the ISR luminos- 
ity at the tp(2S) mass; emc = 0.059 as the detection ef- 
ficiency from simulation with corrections and systematic 
error discussed in Sec. VF, and a 0.90 correction factor 
we obtain 



-DV(2S)^47r27r'' " O^int 



N{tj}{2S) 47r27r0) 

dC/dE ■ CMC 
19.1 ±5.6 ±2.2 nbMeV 
N{ij{2S) -> 47r27r0) • m. 



67r2 • dC/dE ■ euc ■ C 
(1.12 ± 0.33 ± 0.13) X 10^2 keV. 



Using the value of vt^"^^^ from the PDG [23] we calculate 
5v(2S)^4,.27rO = (5.3 ± 1.6 ± 0.6) X lO'^. This is the first 
measurement of that branching ratio. 



VIII. THE J/V PRODUCTION RATE AND 
CONTINUUM 



where cmc = 0.060 is the detection efficiency from 
simulation with the corrections and error discussed in 
Sec. VF. The subscript "471277°" for branching fractions 
refers to the 2(7r+7r~7r°) final state exclusively. 

Using VU"'' = 5.40 ± 0.18 keV [23], wc obtain the 
result = (1.65 ± 0.10 ± 0.18) x 10^^ No 

entry for this branching fraction exists in PDG [23] . As 
was noted in Sec. V G 170 ± 24 events arc observed in 



the J/tp 
value, 



"77° channel. This corresponds to the 



-'^J/i/j— fa;7r+7r~7r*^ * ee 



(2.2 ±0.3 ±0.2) X 10-2 keV 



including the decay rate for uj tt^tt tt . The cor- 
responding fraction is Bjj^^^^+^-^o = (0.41 ± 0.06 ± 
0.04) X 10"^, which is the first measurement of this 
branching fraction. There are also 13 ± 3.6 events from 
the J/ij) ^ LOT] decay. Taking into account the decay 
rates of u and r] to tt+tt-tt^, we obtain Bj/^^^^ = 
(1.47±0.41±0.15) X 10"^, in agreement with the current 



PDG value B 



J/tp—^LOl] 



(1.58 ±0.16) X 10- 



The tp{2S) peak can partly correspond to the decay 
chains V(25) ^ J/ijTT^iT^ 2{tt+tt-tt°) or i/)(2S') ^ 
J/-tpTT+TT~ 2(7r+7r-7r°). Using -Bj/^^2(7r+7r-) = 
(3.61 ± 0.37) X 10"^ from Ref. [8], and assuming that 



B 



+^-^Oj^o is at the same level, as confirmed by 



J/V'— >7r 

a preliminary study of that channel, we estimate that 
43 ± 15 events from the above decay chains contribute 
to the ip{2S) peak. The total number of tl>{2S) signal 
events extracted from a fit to a Gaussian distribution 
for the ip{2S) peak and a polynomial function for the 
continuum is 128 ± 20. The 43 ± 15 events estimated 
above should be subtracted from this total to obtain an 
estimate for the number of direct ipC^S) — > 2(77+77-77°) 
decays. 



Having measured products of a^J^ and branching frac- 
tions for J/ip decaying to 677 and 2K477, it is interesting 
to compare them with the non-resonant cross sections 
(continuum) at that energy. Using a second order poly- 
nomial approximation of the cross sections from Tables I, 
V, and III around the J/tp peak within ±0.4 GeV (events 
from the peak are excluded), the following cross sections 
are obtained for the continuum at the J/ip mass: 

ae^ = 0.54 ±0.03 nb 
O-4,r27r0 = 2.17 ±0.04 nb 
CT2K47T = 0.54 ± 0.04 nb . 



Table VI presents the ratios Bj/^^f ■ afj^ /a^- 
f = 677, 2i4r477. In these ratios most of the experimental 
systematic errors cancel. Also shown are the ratios for 

the J/tp final state [6], and J/tjj 2(77+77-), 

J/i) K+K-n+n', and J/i/' ^ K+R-R+K' [8]. 



for 



TABLE VI: Ratios of the J/tp partial production rates to 
continuum cross sections. The result for /U+/l(- is from Ref. [6]. 
The results for 2(77+77-), K+K'-k+'k- and K+ K' K+ K' are 
from Ref. [8]. 



Final state, / 



j,^^f<yU^/ae+.-^f (MeV) 



3(77+77-) 
2(77+77-77°) 
7^+7^-2(77+77- 
2(77+77-) 

7s:+7s:-77 + 77- 

K+K'K+K- 



106 ± 10 
99.1 ±6.5 

122 ± 10 
85.1 ± 7.9 
166 ± 19 
138 ± 32 
84.12 ±0.67 



The strong decay of the J/tp to an even number of pi- 
ons is forbidden by G-parity conservation and therefore 
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this decay is expected to be dominated by a single pho- 
ton. No such suppression due to G-parity occurs for the 
strong decay of the J /ijj for the other modes. The ratio 
obtained for both of the Btt final states is close to that 
for iX^ n,~ and for 47r, indicating that the single-photon 
exchange dominates for the J /ip decays into these modes. 
For the J/ip decays to final states with kaons, the single- 
photon mechanism may be less dominant, as indicated 
by the significantly larger values of the ratios. 



IX. SUMMARY 

The large 232 fb~^ data sample accumulated near the 
T{AS), together with the excellent resolution, charged 
particle identification, and open trigger of the BABAR 
detector, affords a unique opportunity to study the 

3(7r+7r-), 2(7r+7r-7r°), and K+ K-2{-k+tt-) final states 
produced at low effective e+e~ cm. energy via ISR. 
Not only are the data samples extraordinarily large and 
well measured, but they do not suffer from the relative 
normalization uncertainties which seem to plague direct 
measurements from earlier storage rings such as DCI and 
SPEAR. 

Since the luminosity and efficiencies are understood 
within 3-5% accuracy, these data allow useful measure- 
ments of R, the ratio of the hadronic to di-muon cross 
sections, in the low e^e~ energy regime, providing im- 
portant input for calculating corrections needed for muon 
g — 2 measurements. Cross section measurements for the 
reactions e+e~ — > 3(7r+7r~) and e+e~ — > 2(7r+7r~7r°), 
with about 6% and 10% total systematic errors, respec- 
tively, are significant improvements over earlier measure- 
ments, while no earlier mcasm-ements exist for the reac- 
tion e+e" K+K-2{'K+TT-). 

These final states also provide new information on 
hadron spectroscopy. Structure is observed around 
1.9 GeV/c2, in both the 3(7r+7r-) and 2(77+77-7:°) final 
states, which is similar to that previously observed by the 
DM2 and FOCUS experiments. A simple coherent model 
fits the structure well, and leads to the following param- 
eters: m = 1.87 ±0.02 GeV/c2, T = 0.14 ± 0.03 GeV 
with — 3± 15° relative phase to the continuum. However, 
these values are model dependent. 

Resonance substructure is also seen in the intermediate 
mass combinations for these final states. The tt^tt~ de- 
cay of the p{770) dominates the substructure observed in 
the 3(77+77-) final state, with little evidence for any other 
structures. A model with a single p{770) describes the 
observed distributions very well. On the other hand, not 
only does the 2(77+77-77") final state display contributions 
from p(770)477 but also from ^(782)77+77-77° and w(782)r/ 
final states. Structures that may correspond to /o(980), 
/o(1370) and/or /2(1270) are also visible in intermedi- 
ate mass combinations. A detailed understanding of this 
complex structure will probably require a full PWA using 
all of the six-pion final states. 

Except for a structure around 1.6 GeV, the ratio of 



e+e- 2(77+77-77°) to e+e- — > 3(77+77-) cross sections 
is equal to 3.98 ± 0.06 ± 0.41 over a wide range of cm. 
energies. However, it is no longer flat and substantially 
smaller when the contribution from a;(782) decays to 
77+77-77° is excluded. These results can indicate and help 
to understand more complex structures existing in the 
six-pion final states. 

A resonance-like structure is seen in the e+e- 
a;(782)?7 cross section around 1.6 GeV. A resonance fit 
gives m = 1.645 ± 0.008 GeV/c^, F = 0.114 ± 0.014 GeV 
and peak cross section ctq = 3.08 ± 0.33 nb. This ob- 
ject might be associated with the a)(1650), for which a 
a;(782)?7 final state has been seen [23], but the width 
measured here is substantially smaller. The observed pa- 
rameter values are actually much closer to those of the 
0(1680) [23], but there are no observations for an a;(782)?7 
decay mode of the 4>{1680). 

The energy dependence of the cross section for the 
reaction e+e- K~^K~2{tt'^tt~) from threshold to 
4.5 GeV is measured with about 15% systematic un- 
certainty. There is clear evidence for resonance sub- 
structure including the iir*°(892) and (j). 

These data also allow a study of J/ip and ip{2S) pro- 
duction, and the measurement of the product of decay 
branching fractions and the e+e- width of the J/tjj with 
the best accuracy to date. The results are as follows: 

-Bj/V;^3(7r+7r-) ' ^U"^ 

= (2.37 ± 0.16 ± 0.14) X 10-2 keV , 

Bj/^^K+K-2(tt+tt-) ■^il^ 

= (2.75 ± 0.23 ± 0.17) X 10-2 keV , 
Bj/,p^^(i02n)2(7T+7T-) • B^^x+K- • rf/''' 

= (0.47 ± 0.09 ± 0.03) X 10-^ keV , 

-Bj/^^2(7r+7r-7rO) ' Tg/* 

= (8.9 ± 0.5 ± 1.0) X 10-2 keV , 

J/i/j— >aJ7r+7r~7r^ * ee 

= (2.2 ± 0.3 ± 0.2) X 10-2 keV , 

-B^(2S)^2(7r+7r-7rO) ' T^i^"^^ 

= (1.12 ± 0.33 ± 0.12) X 10-2 keV , 

D Fi/'(2S) 
Otl,(2S)^K+K-2(TT+TT-) ■ J- (le 

= (4.4 ± 2.1 ± 0.3) X 10-3 keV . 

Dominance of the single-photon decay-mechanism for 
J/^f) 3(77+77-) and J/tp — !■ 2(77+77^77°) decays has been 
demonstrated by comparison with the continuum cross 
sections for e+e- — > 3(77+77-) and e+e- 2(77+77-77°). 
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